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Abstract. Shifted surface spline is a frequently used radial function for scattered data interpo- 
lation. The most frequently used error bounds for this radial function are the one raised by Wu and 
Schaback in [T7] and the one raised by Madych and Nelson in [TJ. Both are 0(d l ) as d — > 0, where 
I is a positive integer and d is the well-known fill-distance which roughly speaking measures the 
spacing of the data points. Recently Luh raised an exponential- type error bound with convergence 
rate O(w^) as d — > where < u < 1 is a fixed constant which can be accurately computed. 
Although the exponential-type error bound converges much faster than the polynomial-type error 
bound, the constant u> is intensely influenced by the dimension n in the sense to — > 1 rapidly as 
n — > oo. Here the variable x of both the interpolated and interpolating functions lies in R n . In this 
paper we present an error bound which is 0(^fdu)'&) where < u' < 1 is a fixed constant for any 
fixed n, and is only mildly influenced by n. In other words, u)' — ► 1 very slowly as n — + oo, and 
lu' << oj, especially for high dimensions. Moreover, w' can be accurately computed without slight 
difficulty. 
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1 Introduction 



In the theory of radial basis functions, it's well known that any conditionally positive definite radial 
function can form an interpolant for any set of scattered data. We make a simple sketch of this 
process as follows. 

Suppose h is a continuous function on R n which is strictly conditionally positive definite of order 
m. For any set of data points (xj, fj),j = 1,...,N, where X — {x\, . . . , xn} is a subset of R n and 
the fjS are real or complex numbers, there is a unique function of the form 

JV 

s(x) = p(x) + Cjh(x — xj) (1) 
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where p{x) is a polynomial in P!^ l _ l , satisfying 



N 

5^c,-g(x J )=0 (2) 



for all polynomials g in P^-i an( i 



A' 



+ ^ c,-ft(a;t - Xj) = fi, i = 1, . . . ,N (3) 

if X is a determining set for Pj£-i- 

A complete treatment of this topic can be seen in |13j and many other papers. 

The function s(x) is called the /i-spline interpolant of the data points and is of central importance 
in the theory of radial basis functions. In this paper h always denotes a radial function in the sense 
that the value of h(x) is completely determined by the norm \x\ of x. Here, P^-i denotes the class 
of those n- variable polynomials of degree not more than m — 1. 

In this paper we are mainly interested in a radial function called shifted surface spline defined 

by 

h(x) := (-l) m (M 2 + c 2 )?log(\x\ 2 + c 2 )^, A e Z+, m = 1 + ^, c> 0, 

x e R n , A, n even, (4) 

where \x\ is the Euclidean norm of x, and A, c are constants. In fact, the definition of shifted surface 
spline covers odd dimensions. For odd dimensions, it's of the form 

h(x) := {-l)^ x -^(\x\ 2 +c 2 ) x -^, nodd, AeZ+ = {1,2,3,...} 

TL 

and A > — . (5) 

2 v ; 

However, this is just multiquadric and we treat it in another paper |12j . Therefore we will not discuss 
it. Instead, we will focus on (4) and even dimensions only. 



1.1 Polynomials and Simplices 

Let E denote an n-dimensional simplex]!] with vertices V\, . . . , u n +i. If we adopt barycentric coor- 
dinates, then any point ie£ can be written as a convex combination of the vertices: 

n+l n+ 1 

X = ^2 ^i v ii J! ^ = ^ ^ - " 
i=l i=l 

We define the "equally spaced" points of degree k to be those points whose barycentric coordinates 
are of the form 

(k\/k, fe/fc, . . . , k n+ i/k), ki nonnegative integers and ki + • • • + k n +\ = k. 

It's easily seen that the number of such points is exactly dimPJ} , i.e., the dimension of PJ} . In this 
section we use N to denote dimP£. 

The above-defined equally spaced points can induce a polynomial interpolation process as follows. 
Let Xi, . . . , Xn be the equally spaced points in E of degree k. The associated Lagrange polynomials 
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li of degree k are defined by the condition h(xj) — 8ij, 1 < i, j < N . For any continuous map 
/ € C(E), (Ukf)(x) := J2iLi f(xi)U(x) is its interpolating polynomial. If both spaces are equipped 
with the supremum norm, the mapping 

II, : C(E) - P£ 

has a well-known norm 

N 

\\U k \\ = max^2\h(x)\ 

i=l 

which is the maximum value of the famous Lebesgue function. It's easily seen that for any p £ PJ} , 

HpIIoc := max|p(x)| < ||IT fe || max \p(xi)\. 

xGE 1<i<N 

The next result is important in our construction of the error bound, and we cite it directly from [SJ. 

/ 2k — 1 \ 

Lemma 1.1 For the above equally spaced points {x\, . . . , Xn}, \\^-k\\ ^ ( ^ )• Moreover, as 
n->oo, 1 1 IXjfe 1 1 -> f 2fc fc 1 V 

Then we need another lemma which must be proven because it plays a crucial role in our development. 

Lemma 1.2 Let Q C i?" 6e aw n simplex in R n and Y be the set of equally spaced points of degree 
k in Q. Then, for any point x in Q, there is a measure a supported on Y such that 

J P(y)dcr(y) = p(x) 

for all p in PJ} , and 

V)- 

Proof. Let Y = {yi, . . . ,yN} be the set of equally spaced points of degree k in Q. Denote PJ? 
by V. For any x £ Q, let S x be the point-evaluation functional. Define T : V — > T(V) C R N 
by T(v) = (6 Vi (v)) yieY - Then T is injective. Define ip on T(V) by ip(w) = S x (T- 1 w). By the 
Hahn-Banach theorem, ip has a norm-preserving extension i[> ex t to R . By the Riesz representation 
theorem, each linear functional on R N can be represented by the inner product with a fixed vector. 
Thus, there exists z S R N with 

N 

4>ext(w) = ^ z i w i 

J'=l 

and ||2||(ijjv)* = HV'extll- If we adopt the ioo-norm on R N , the dual norm will be the Zi-norm. Thus 
\\z\\ {R N ) y = \\z\\ 1 = U ext \\ = U\\=\\5 x T-l\\. 

Now, for any p £ V, by setting w = T(p), we have 

N N 

S x (p) = 5 x (T^ 1 w) = ^(w) = i) ext (w) = ^ z j w j = X! z 3 S vAp)- 

3=1 0=1 
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This gives 



where \z x \ + ■■■ + \z N \ = \\S x T- l \ 
Note that 



Therefore \z\\ H h \z N \ < 



< 



< 



N 



3=1 



< 

2k- 1 
k 



sup 
w e T(V) 

sup 

w e t(v) 

to/0 
sup 



|^r-»H 
\s x p\ 

\ t (p)\\r n 
\p(x)\ 



max j= i,..., w |p(%)| 



n fe || max 1= 



■j = h-,N 



\p(y 3 )\ 



(6) 



sup - 

2k- 1 
k 



and our lemma follows immediately from (6) by letting 



1.2 Radial Functions and Borel Measures 

Our theory is based on a fundamental fact that any continuous conditionally positive definite radial 
function corresponds to a unique positive Borel measure. Before discussing this property in detail, 
we first clarify some symbols and definitions. In this paper V denotes the space of all compactly 
supported and infinitely differentiable complex- valued functions on R n . For each function </> in T>, 
its Fourier transform is 

to -/«-«■ 

Then we have the following lemma which is introduced in [8] but modified by Madych and Nelson 
in [H]. 

Lemma 1.3 For any continuous conditionally positive definite function h on R n of order m, there 
are a unique positive Borel measure \x on R n ~ {0} and constants a r , \r\ — 2m such that for all 



h(x)ip(x)dx 



\r\<2m 
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+ E ( 7 ) 

|r| <2m 

, where for every choice of complex numbers c a , \a\ — m, 

Here x is a function in T> such that 1 — x(£) has a zero of order 2m + 1 at £ = 0; &o£/i 0/ the integrals 

0<\(\<1 J\(\>1 

are finite. The choice of x affects the value of the coefficients a r for \r\ < 2m. 

2 Main Result 

In order to show our main result, we need some lemmas, including the famous Stirling's formula. 
Stirling's Formula: nl ~ \ / 2nn(j) n . 

The approximation is very reliable even for small n. For example, when n = 10, the relative 
error is only 0.83%. The larger n is, the better the approximation is. For further details, we refer 
the reader to [5] and [6]. 

Lemma 2.1 For any positive integer k, 



vW <2 k 
k\ - • 

Proof. This inequality holds for k = 1 obviously. We proceed by induction. 



y/[2(fc + l)]! y/(2fc + 2)! _ y/(2k)\ v /(2fc + 2)(2fc + l) 

(k + l)\ ~ k\{k + l) k\ k + 1 

< TP)! V(2A: + 2)2 fc (2k + 2) _ 

k\ k + 1 ~ k + 1 ' H 

Now recall that the function h defined in (4) is conditionally positive definite of order m = 1 + -| . 
This can be found in [3] and many relevant papers. Its Fourier transform [7] is 

h(6) = l{X,n)\e\- x - n iC^{c\9\) (8) 

where Z(A, n) > is a constant depending on A and n, and JC v (t) — t^JC y (t), K, v (t) being the modified 
Bcssel function of the second kindp]. Then we have the following lemma. 

Lemma 2.2 Let h be as in (4) and m be its order of conditional positive definiteness. There exists 
a positive constant p such that 

< l(\, n) ■ y| • n ■ a n ■ c x ~ k ■ A • p k ■ k\ (9) 

for all integer k > 2m + 2 where fi is defined in (6), a n denotes the volume of the unit ball in R n , c 
is as in (4), and Aq is a positive constant. 
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Proof. We first transform the integral of the left-hand side of the inequality into a simpler form. 

\z\ k i(\, n)£»±x(cKi)Kr A - n de 



l(X,n)c^ / \^\ k -^ -JC^(c\^ 



j(A, n)c^ . / j / |C| fc -^ • * - dg 



1+A /7T / ju »+a r 



n-l 



= i(A,n)c > ■^■"■*nJ o r * -^T^dr 
= i(A,n)c > ^-.n-a„.-^y o <*r 

= l(X,n)c * ^-.n.^.-^.^^J^ ~^—dr 

x \-k f°° rk \ , ,/ , n-A-3 

= /(A, n)\ — ■ n ■ a n ■ c / — ar where k =k + 
V 'V 2 " Jo e r 



2 

Note that fc>2m + 2 = 4 + A implies fc' > ™+|+5 > . 

Now we divide the proof into three cases. Let k" — \k'~\ which is the smallest integer greater 
than or equal to k! . 

Casel. Assume k" > k. Let k" = k + s. Then 

roo k' roo k" 

J —dr<J —dr = k"\ = (k + s)(k + s-l)---(k + l)k\ 

and 

roo fe'+l roc k" + 1 

/ —--dr< — dr = (k" + 1)1 = (k + s + l)(k + s) ■■■(k + 2)(k + l)kl. 

Jo eT Jo eT 

Note that 

{k + s + l)(k + s) ■ ■ ■ (k + 2) _ k + s + 1 
(fc + s)(fc + s - 1) • • • (fc + 1) ~~ fe + 1 
The condition k > 2m + 2 implies that 

k + s + 1 2m + 3 + s 1 s 

< = H . 

k + l ~ 2m + 3 2m + 3 

/•oo k" + l 

/ < A -p k+1 ■ (k+ l)! 

if J* °° < A • p k ■ hi. The smallest k" is k' ' = 2m + 2 + s when k = 2m + 2. Now, 

f°° r k « 

/ dr = k' Q '\ = (2m + 2 + s)(2m + l + s) ■ ■■ (2m + 3) (2m + 2)! 

Jo e r 
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(2m + 2 + S )(2m + l + S )...(2m + 3) 2m+2 

= A -p 2m+2 -(2m + 2)! 

(2m + 2 + s)(2m + 1 + s) • • • (2m + 3) 



where A = 



p2m+2 



It follows that J °° ^-dr < A • p fc • fc! for all fc > 2m + 2. 

Case2. Assume fc" < fc. Let fc" = k — s where s > 0. Then 



/■oo fc' roo fc" 

/ — dr< / rf r = fc"! = (fc-s)! = — - r • fc! 



fc(fc-l)---(fc-s + l) 
and 



/ -^dr < / 



-dr 



e' 



Note that 



«=" + = (* - + D' = (t + 1)t .V» + 2) ■ * + ')'• 



Let p = 1. Then 



(fc + l)fc • • • (fc - s + 2) ' fc(fc - 1) • • • (fc - s + 1) 

fc(fc-l)---(fc-5+l) 

(fc + l)fc---(fc-s + 2) 
(k-s + l) 
fc + 1 

< 1. 



/■oo fc"+l 

/ dr < A -p fe+1 • (fc+ 1)! 

Jo e r 



if J °° ^—^dr < A a ■ p k ■ fc!. The smallest fc is fco = 2m + 2. Hence the smallest fc" is fc ' = fc — s = 
2m + 2 - s. Now, 



/ — dr = fc^! = (2m + 2-s)! = (fc -s)! 
Jo e 



! (fco!) 



fc (fc -l)---(fco-s + l) 
= A • p k ° ■ fc ! where A = 



(2m + 2) (2m + !)■■■ (2m - s + 3) 



It follows that J °° ^rdr < A • p k ■ fc! for all fc > 2m + 2. 
Case3. Assume fc" = fc. Then 

roo fc' fOO fc" /■OO fe' + l 



/■OO fc' /-OO fc" /-OO fc'+l 

/ — dr< / dr = kl and / dr<(fc + l)!. 

Jo e r Jo e r J e r 
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Let p = l. Then f °° ^-rdr < A • p k ■ fc! for all fc where A = 1. 

The lemma is now an immediate result of the three cases. jj 

Remark: For the convenience of the reader we should express the constants Ao and p in a clear 
form. It's easily shown that 

(a) fc" > fc if and only if n — A > 3, 

(b) fc" < fc if and only if n — A < 1, and 

(c) fc" = fc if and only if 1 < n — A < 3, 

where fc" and fc are as in the proof of the lemma. We thus have the following situations. 

(a) n - A > 3. Let s = \ n ~^~ 3 ]. Then 

s , . (2m + 2 + s)(2m+l + s)---(2m + 3) 
p = H and A = - - — - — - - -. 

(b) n - A < 1. Let s = - f^^l^l • Then 



p = 1 and Aq 



1 



{2m + 2) (2m + 1) ■ • • (2m - s + 3) ' 
(c)l < n- A < 3. We have 

p = 1 anrf Ao = 1. 

Before introducing our main theorem, we must introduce a function space called native space, 
denoted by Ch.mj for each conditionally positive definite radial function h of order m. If 

V m = {(f> G D : y x a (j)(x)dx = /or ct/Z |a| < m} 

, then C/j. m is the class of those continuous functions / which satisfy 



f (x)4>(x)dx 



< c(f) h(x- y)cj)(x)cj)(y)dxdy (10) 



1/2 



for some constant c(/) and all <j> in T> m . If / € Cft, jm , let denote the smallest constant c(/) for 
which (10) is true. Then || • ||^ is a semi-norm and Ch. m is a semi-Hilbert space; in the case m = it 
is a norm and a Hilbert space respectively. For further details, we refer the reader to [T5] and [Ti] , 
This definition of native space is introduced by Madych and Nelson, and characterized by Luh in 
[S] and [TU]. Although there is an equivalent definition [IB] which is easier to handle, we still adopt 
Madych and Nelson's definition to show the author's respect for them. 
Now we have come to the main theorem of this paper. 

Theorem 2.3 Let h be as in (4), and bo be any positive number. Let be any subset of R n satisfying 
the property that for any x in f2, there is an n simplex E of diameter equal to bo and x G E C Q. 
There are positive constants <5o,ci,a/, < a/ < 1, for which the following is true:Lf f G Ch.m, the 
native space induced by h, and s is the h spline that interpolates f on a subset X of R n , then 

\f{x)-s{x)\ < Cl Vj(a/)*- \\f\\ h (11) 

for all x in fl and < 5 < So if & satisfies the property that for any x in Q and any number r with 
3^7 < t < bo, there is an n simplex Q with diameter diamQ=r. x G Q C tt, such that for any integer 
k with -^jfi < fc < ^f, there is in Q an equally spaced set of centers from X of degree k-l.(Once 8 
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and k are chosen, Q is in fact an n simplex of diameter k5 and X can be chosen to consist only of 
the equally spaced points in Q of degree k-1.) The number C is defined by 



C:=maxjV, — j, pf = Z 

where p and c appear in Lemma2.2 and (4) respectively. Here \\f\\h is the h-norm of f in the 
native space. The numbers Sq,ci and uu' are given by So := 3^,^+1) w ^ ere m appears in (4); 

c\ := y/l(X, n) • (tt/2) 1 / 4 • y/na n ■ c A / 2 • \J Ao\/3C • \J (167r) _1 where A is as in (4), l(X,n) appears in 
(8), a n is the volume of the unit ball in R n , and Ao, together with the computation of p, is defined 
in Lemma2.2 and the remark following its proof; oj' := (^) 1 ^ 3C ■ 

Proof. Let So, and C be as in the statement of the theorem. For any < 6 < So, we have 
< ^ — 3C(tk+i) anc ^ < 3C<5 < ^zfi- Since < 1, there exists an integer k such that 

1 < 3C*<Sfc < 2. 

For such k, Sk < ^ <b Q , ^ < k < !f, and 8p'Sk < §. 

For any x £ Q, choose arbitrarily an n simplex Q of diameter diamQ = Sk with vertices 
vq, vi, . . . , v n such that x £ Q C O. Let xx,...,xn be equally spaced points of degree k-1 on 
Q where N = dimPf}^. By (9) and Lemma2.1, whenever k > m, 

< 0(A^)-(^/2) V4 - V^-c A/2 -c" fe - SK~ -{2p) k . (12) 
Recall that Theorem4.2 of [H] implies that 

\f(x) - s(x)\ < c k \\f\\ H ■ [ \y - x\ k d\a\(y) (13) 

JR" 

whenever k > m, and a is any measure supported on X such that 

p(y)do-(y) = p(x) (14) 



for all polynomials p in P£_ 1 - 

Let a be the measure supported on {x\, . . . , xn} as mentioned in Lemmal.2. We essentially need 
to bound the quantity 

I =c k \y- x\ k d\a\(y) 



only. Thus for k mentioned as above and Ao defined in Lemma2.2, by Lemmal.2, 

/ < ^l^){n/2) 1 ' i ^ l c x ' 2 c- k ^o{2p)\5k) k (^ 2{k k ^Y 1 ^ 

~ V / Z(A,n)(7r/2) 1/4 V^Q~c A/2 c~ fc v /A^(2p) fc ((5fc) fc ^ - 1 2 2(fc ~ 1} by Stirling's Formula 

V 71 " yk — 1 

~ ^l(X,n)(n/2) 1/i ^^c x/2 c- k ^/Ao'(2p) k (Sk) k ^=^4 k - 1 when k is large 

V 71 " vfc 
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(tt/2) 1 /4 aA ^ c a/2 v ^. 
( 7 r/2) 1/4 % A^c A / 2 VA^- 
(tt/2) 1 / 4 V^c a / 2 VA^- 
(7r/2) 1/4 V^c A / 2 v/A^- 



167T 
1 I 

/ 167T VA? V 3 , 
1 1 /2\^ 



/ I6tt y/k V 3 
1 1 



; 167r \fk 
1 1 , 



2\ 3C 
3 



— = |q/l 7 where u)' = ( — 
^n/I V 3 - 



(tt/2) 1/4 V^c A/2 v^-^VsCJ [a/] * 

V 16"7T 



( 7 r/2) 1 / 4 V^c A/2 



T6tT 



Our theorem thus follows from (13). ft 

Remark. In the preceding theorem we didn't mention the well-known fill-distance. In fact 5 is in 
spirit equivalent to the fill-distance d{Q, X) :— sup yg Q mini<i<7v \\y — where X — {xi, . . . , xn} 
as mentioned in the proof. Note that S — > if and only if d(Q,X) — > 0. However we avoid using 
fill-distance because in our approach the data points are not purely scattered. This to some ex- 
tent seems to be a drawback. However it does not pose any trouble for us both theoretically and 
practically. The equally spaced centers 27, . . . , xjq in the simplex Q are friendly and easily tractable. 



3 Comparison 

The exponential-type error bound for (4) is presented by Luh in [11] and is of the form 

|/(a:)-a(a;)|<ciw*||/|| fc (15) 



where c x = y/l(X, n)(ir /2) 1 / i y / na„c A ' 2 \/Ao, S is equivalent to the fill-distance and 



3, 



where 



C = max|2p'V^e 2 ™ 7 ", 



' 3& J ' 



, p and c being the same as this paper, bo be the side length of a cube, and j n being defined 
recursively by 

7i = 2, 7„ = 2n(l + 7„-i) if n > 1. 
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The constant c\ is almost the same as the c\ in (11). The number 60 plays the same role as the 60 
of Theorem2.3. However, 7 „ — > 00 rapidly as n — > 00. This can be seen by 

7l = 2, 72 = 12, 73 = 78, 74 = 632, 7s = 6330, • • • 

The fast growth of 7 „ forces e 2 ™ 7 ™ and hence C to grow rapidly as dimension n — > 00. This means 
that the crucial constant in (15) turns to 1 rapidly as n —* 00, making the error bound (15) 
meaningless for high dimensions. 

The advantages of our new approach are:first, there is \/S in (11) which contributes to the 
convergence rate of the error bound as 5 — > 0; second, the crucial constant u>' in (11) are only mildly 
dependent of dimension n. Although u>' dependends on p which in turn depends on n, the situation 
is much better. In fact, u>' can be made completely independent of n by changing A in (4) to keep 
n — A < 3. This can be seen in the remark following Lemma2.2. In other words, we have significantly 
improved the error bound (15). 

References 

[1] M. Abramowitz and I. Stegun, A Handbook of Mathematical Functions, Dover Publications, 
New York, (1970). 

[2] L.P. Bos, Bounding the Lebesgue Function for Lagrange Interpolation in a Simplex, J. of Approx. 
Theory 38, (1983), 43-59. 

[3] N. Dyn, Interpolation and Approximation by Radial and Related Functions, Approximation 
Theory VI,(C.K. Chui, L.L. Schumaker and J. Ward eds.), Academic Press, (1989), 211-234. 

[4] W. Fleming, Functions of Several Variables, Second Edition, Springer- Verlag, (1977). 

[5] J. von zur Gathen and J. Gerhard, Modern Computer Algebra, Second Edition, (2003), Cam- 
bridge University Press. 

[6] R.L. Graham, D.E. Knuth, and O. Patashnik, Concrete Mathematics, Addison- Wesley, Reading 
MA, 2nd cdition(1994). 

[7] I.M. Gelfand and G.E. Shilov, Generalized Functions, Vol.1, Academic Press, (1964). 

[8] I.M. Gelfand and N.Y. Vilcnkin, Generalized Functions, Vol.4, Academics Press, New York, 
(1964). 

[9] L.T. Luh, The Equivalence Theory of Native Spaces, Approx. Theory Appl. (2001), 17:1, 76-96. 

[10] L.T. Luh, The Embedding Theory of Native Spaces, Approx. Theory Appl. (2001), 17:4, 90-104. 

[11] L.T. Luh, The High-level Error Bound for Shifted Surface Spline Interpolation, in review. 

[12] L.T. Luh, An Improved Error Bound for Multiquadric Interpolation, in review. 

[13] W.R. Madych and S.A. Nelson, Multivariate interpolation and conditionally positive definite 
function, Approx. Theory Appl. 4, No. 4(1988), 77-89. 

[14] W.R. Madych and S.A. Nelson, Multivariate interpolation and conditionally positive definite 
function, II, Math. Comp. 54(1990), 211-230. 



11 



[15] W.R. Madych and S.A. Nelson, Bounds on Multivariate Polynomials and Exponential Error 
Estimates for Multivariate Interpolation, J. Approx. Theory 70, (1992), 94-114. 

[16] H. Wendland, Scattered Data Approximation, Cambridge University Press, (2005). 

[17] Z. Wu and R. Schaback, Local Error Estimates for Radial Basis Function Interpolation of 
Scattered Data, IMA J. of Numerical Analysis, (1993). 



12 



